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Subpoint  Prediction  for  Direct  Readout 
Meteorological  Satellites 

L.  E.  EBER' 


ABSTRACT 

The  National  Environmental  Satellite  Service  (NESS)  provides  orbital  information  on  me- 
teoroloKical  satellites  with  direct  transmission  systems,  through  APT  (Automatic  Picture  Trans- 
mission) Predict  messages  sent  over  standard  weather  communications  networks.  With  peri- 
odic access  to  this  information,  operators  of  independent  APT  ground  receiving  stations  can 
extrapolate,  by  means  of  nodal  period  and  nodal  increment,  to  determine  future  orbits  within 
receiving  range  of  their  station.  A  technique  for  the  prediction  of  subpoint  location  along  an 
orbit  as  a  function  of  time  after  ascending  node  was  developed  from  consideration  of  Kepler's 
laws  and  derived  expressions  for  the  force  due  to  the  earth's  gravitational  potential.  Subpoint 
latitudes  and  longitudes  computed  by  this  technique  are  within  0.1  degree  of  those  given  in  NESS 
predictions. 


INTRODUCTION 

Users  of  APT  (Automatic  Picture  Transmis- 
sion) and  Direct  Readout  Scanning  Radiometer 
data  from  meteoroio,<?ical  satellites,  who  operate 
APT  ground  receiving  stations,  need  orbital  in- 
formation for  scheduling  transmission  pickup  and 
for  tracking  and  locating  subpoints  for  data  iden- 
tification. Such  information  is  provided  by  the 
National  Environmental  Satellite  Service  (NESS) 
through  APT  Predict  messages  transmitted  over 
meteorological  teletype  circuits.'  In  cases  where 
these  messages  are  not  regularly  available,  or 
where  advance  scheduling  is  desired,  the  dates, 
times,  and  longitudes  of  the  ascending  nodes  for 
selected  orbits  can  be  extrapolated  by  means  of 
nodal  period  and  nodal  increment  as  much  as  a 
month  ahead  with  acceptable  accuracy.  Subpoint 
locations  at  selected  time  invervals  along  the 
track  also  can  be  predicted  for  future  orbits  if 
certain  other  orbital  characteristics  are  known. 


'  Southwest  Fislieries  Center,  National-  Marine  Fislieries 
Service.  NOAA,  La  Jolla,  CA  92037. 

-ESSA  (Environmental  Science  Services  Administration) 
Direct  Transmission  System  Users  Guide,  prepared  by  the 
National  Environmental  Satellite  Center,  U.S.  Department 
of  ('(immerce,  19()9. 


These  include  the  semimajor  axis  and  eccentric- 
ity of  the  orbital  ellipse,  the  orbital  inclination, 
and  the  geocentric  angle  of  the  perigee  measured, 
at  a  known  reference  time,  from  the  ascending 
node  in  the  orbital  plane. 


UNDISTURBED  ORBITAL  MOTION 

The  orbit  of  a  meteorological  satellite  is  an  el- 
lipse with  the  earth's  center  at  one  focus.  Accord- 
ing to  Kepler's  laws,  the  motion  of  a  satellite  is 
such  that  a  ray  from  the  focus  to  the  satellite 
sweeps  out  equal  areas  in  equal  times.  Conse- 
quently the  angular  velocity  of  the  satellite  is 
greatest  at  perigee  when  the  satellite  is  at  its 
nearest  approach  to  the  earth  and  least  at  apogee 
when  it  is  farthest  away.  The  area  swept  out  per 
unit  time  by  the  satellite  can  be  expressed  by: 


i?2rf 


TiAB 


2dt 


where  d  a/dt  is  the  angular  velocity,  R  is  the  dis- 
tance from  the  center  of  the  earth  to  the  satel- 
lite, A  and  B  are,  respectively,  the  semimajor 
and  semiminor  axes  of  the  orbital  ellipse,  and  P 
is  the  orbital  period.  Rearrangement  of  terms 


gives  an  explicit  equation  for  the  angular  veloc- 
ity, as  follows: 


d  a 
dt 


2uAB 


(1) 


The  values  of  A  and  B  are  related  to  the  eccen- 
tricity, e,  of  an  ellipse  by  the  expression 


_F_ 
A 


where 


F  =yJA^   -52 


and  represents  the  distance  from  the  center  to 
either  focus  of  the  ellipse. 

The  relation  between  R  and  a  can  be  obtained 
from  the  basic  equation  of  an  ellipse  in  polar  co- 
ordinates, with  the  origin  at  one  focus  and  the 
polar  axis  coincident  with  the  major  axis  of  the 
ellipse: 


R  = 


A   (1 


1  ±  e  cos  a 


(2) 


DISTURBED  ORBITAL  MOTION 

Because  the  earth  is  not  a  concentrically  homo- 
geneous sphere,  but  rather  an  oblate  ellipsoid 
of  revolution,  the  angular  velocity  is  not  exactly 
defined  by  Equation  (1).  The  elliptic  motion  is 
governed  by  a  force  which  can  be  defined  in 
terms  of  the  Newtonian  gravitational  potential 
U,  as  follows  (Menzel,  1960): 


U  = 


IG  M      K^  (I  -J)  (I  -3sin^  4>) 


R 


+ 


2R^ 


(3) 


where  M  is  the  mass  of  the  earth, 

K^  is  the  constant  of  gravitation, 

I  is  the  moment  of  inertia  about  the  polar 

axis, 
J  is  the  moment  of  inertia  about  any  dia- 
meter in  the  earth's  equatorial  plane, 
and 
4)  is  the  angle  between  the  equatorial 
plane  and  a  line  from  the  earth's  cen- 
ter to  the  satellite  (i.e.,  cf)  =  latitude). 

The  expression  (I  —J)  can  be  related  to  the 
earth's  mass  M,  the  equatorial  radius  a,  the 


earth's  flattening/,  and  the  ratio  q  of  the  earth's 
centripetal  acceleration  at  the  equator  to  gravity, 
as  follows: 


I  -J  =(2/3)  Ma^  if  -q/2). 


(4) 


Substituting  Equation   (4)   in  Equation  (3),   we 
get 


U  = 


K^  M 


R 


+ 


IP  Ma^  {f  -q/2)  (1^3  sin^  (/>) 
3R' 


(5) 


A  more  precise  form  of  expression  for  the 
earth's  gravitational  potential  is  given  in  terms 
of  spherical  harmonics  by 


U  = 


IP  M       IP  M 


R 


R 


QO 

2  J 


n 


-^  I   Pn   (sin  </)) 


where  the  J„  are  constants  which  can  be  deter- 
mined from  observations  of  earth  satellites,  and 
the  Pn  (sin  c^)  are  Legendre  polynomials  (Run- 
corn, 1967).  The  first  term  in  the  series  (corre- 
sponding to  )i  =2)  represents  the  main  effect  of 
the  earth's  oblateness.  The  value  of  J,  has  been 
determined  as  1.0826  XIO'  ,  which  is  400  times 
greater  than  J  for  any  of  the  higher  n.  The 
Legendre  polynomial  for  w  =2  is 

P,  (sin  4>)  =  (3/2)  sin^  </>  - 1/2. 

If  terms  in  higher  n  are  omitted,  the  expression 
for  gravitational  potential  becomes 


U 


IP  M       IP  M 


a2  J,    (1  -3sin2  <J)) 


R 


2R' 


This  expression  is  equivalent  to  Equation  (5) 
if  terms  of  second  and  higher  order  are  omitted 
in  the  following  relation  between  J.  and  the  con- 
stants /  and  q: 

J,  =  (2/3)  (/  -  q/2) 

-(1/3)/ 2    +(l/2)q2    +{2/2l)fq. 


The  values  of  /  and  q  have  been  determined  as 
3.353  X  10  3  and  3.468  X  10  ^  ,  respectively, 
so  that 

(2/3)  ( /  -  q/2)  =  1.0794   X  10  ^  ^  _ 

The  value  of  A'2 Mis  3.986   X  10^   km=  sec"  2. 

The  gravitational  acceleration  is  obtained  by 
differentiating  (5)  with  respect  to  R  to  get 


dU 
dR 


K^  M 
R^ 


K^  Ma^  {f  -q/2)  (1-3  sin^  (^) 
R' 


•(6) 


This  expression  for  gravitational  acceleration  due 
to  the  mass  of  the  earth  could  be  equated  directly 
to  the  centripetal  acceleration  of  an  earth  satel- 
lite if  it  were  traveling  in  a  true  circular  orbit. 
In  an  elliptical  orbit,  however,  there  is  a  devia- 
tion from  circular  motion  which  can  be  identified 
by  substituting  into  (1)  the  following  expression 
from  Kepler's  third  law  for  the  orbital  period: 


2tiA 


3/2 


K  y/M 


The  substitution  gives 


(7) 


from  which  we  obtain  the  following  equation  for 
centripetal  acceleration: 


earth's  gravitational  acceleration  to  the  centri- 
petal acceleration  of  an  earth  satellite  as  follows: 


K(<!^%^(jm 


+ 


/       RA\   R' 

K^Ma^  (f  -q/2)  (1-3  sm'-4>)\ 
R'  J 


and  the  angular  motion  is 

d_a_B_  /  K^M 
dt   ^  R    \  R^A 


K^M  aM  /  -  q/2)  (1-3  sin^  (/>)> 


SUBPOINT  LOCATION 

Satellite  orbit  predictions  are  given  in  terms 
of  subpoint  locations  at  fixed  intervals  of  time 
following  the  northward  equator  crossing,  re- 
ferred to  as  the  ascending  node.  Numerical  inte- 
gration of  Equation  (9),  carried  out  with  suitably 
small  time  increments,  can  be  made  to  yield  suc- 
cessive values  of  a  requisite  for  subpoint  compu- 
tation. If  a  is  taken  to  be  zero  when  the  satellite 
is  at  perigee,  then  the  plus  sign  preceding  e 
cos  a  in  Equation  (2)  is  applicable  and  R  is  given 
by 


R 


A  (1  -e') 
1    +e  cos  a 


(10) 


R 


(I  a 
dt 


IP  M  B^ 
R^  A 


(8) 


The  right  side  of  Equation  (8)  can  be  considered 
as  the  product  of  {B^/RA)  and  the  quantity 
(K^M/R-).  The  latter  is  equivalent,  except  in 
sign,  to  the  first  term  on  the  right  side  of  Equa- 
tion (6)  and  represents  the  gravitational  accel- 
eration for  undisturbed  motion.  The  correction 
for  the  ellipticity  of  an  orbit  is  contained  in  the 
factor  B^/RA.  Its  value  approaches  1  when  an 
orbit  is  nearly  circular.  By  applying  this  factor  to 
the  right  side  of  Equation  (6)  we  can  relate  the 


R  can  be  computed  from  Equation  (10)  for  each 
iteration  in  the  numerical  integration  of  Equa- 
tion (9)  using  the  value  of  a  from  the  preceding 
iteration  increased  by  one-half  the  average 
angular  displacement  [360  (  A  t)/2P\  during  the 
selected  time  increment,  A^-  In  order  to  start 
the  numerical  integration  of  Equation  (9)  at  the 
ascending  node,  an  initial  value  of  a  must  be  de- 
termined from  the  orientation  of  the  orbital  el- 
lipse. The  latter  is  expressed  as  the  argument 
of  perigee,  which  is  the  geocentric  angle  of  the 
perigee  measured  in  the  orbital  plane  from  the 
ascending  node  in  the  direction  of  motion.  The 


Figure  1.  — Right  spherical  triangle  formed  by  the  equator, 
the  orbital  subpoint  track  and  the  meridian  through  the  sub- 
point  ((/),  A). 


initial  value  of  a  is  obtained  by  subtracting  the 
argument  of  perigee  from  360  . 

The  gravity  effect  of  the  earth's  equatorial 
bulge  causes  the  perigee  to  shift  slightly  during 
each  orbit.  Consequently,  it  is  necessary  to  know 
the  rate  at  which  the  orbital  ellipse  is  shifting  in 
order  to  determine  the  correct  initial  values  of 
a  for  prediction  of  future  orbits. 

The  locus  of  orbital  subpoints  for  an  earth 
satellite  traces  a  great  circle  on  a  sphere  concen- 
tric with  the  earth.  Each  subpoint  can  be  repre- 
sented by  a  pair  of  coordinates  (4>,  A),  where  4> 
is  the  geocentric  latitude  and  A  is  the  longitudinal 
displacement  of  the  subpoint  measured  from  the 
longitude  of  the  ascending  node.  Subpoint  co- 
ordinates corresponding  to  any  point  of  the  orbit 
can  be  determined  from  its  angular  displacement, 
measured  along  the  orbital  track  from  the  as- 
cending node.  If  a  is  the  angular  coordinate 
(measured  from  perigee,  as  defined  in  the  pre- 
ceding section)  and  a,,  is  the  value  of  a  at  the  as- 
cending node,  the  angular  displacement  is 
(  «—«„).  Figure  1  shows  these  quantities  as  com- 
ponents of  a  right  spherical  triangle  whose  ver- 
tices are  defined  by  the  subpoint  ((/>,  A),  the  co- 
ordinates of  the  ascending  node  (0,  0)  and  the 
point  (0,  A).  The  two  sides  forming  the  right  angle 
are  (/)  and  A  respectively.  The  opposite  side  is 
(  a  —a„)  and  is  a  segment  of  the  subpoint  track. 
The  angle  opposite  </>,  denoted  by  /,  is  the  inclina- 


tion of  the  orbital  plane  to  the  equatorial  plane. 
Application  of  the  trigonometric  formulas  de- 
rived from  Napier's  rules  yields  the  following  re- 
lationships: 


sin  </)  =  sin(/  )sin(  a  - 

-a,,), 

(11) 

tan  A  =  cos(('  )tan(  a 

^a„],  and 

(12) 

tan  (}>  =  tan(i  )sin  A. 

(13) 

Equations  (11)  and  (12)  give  explicit  formulas  for 
(}>  and  A  as  a  function  of  (/)  and  (  a  — a„).  Equa- 
tion (13)  represents  the  great  circle  traced  by  the 
orbital  subpoints. 

ORBIT  PREDICTION 

The  concepts  developed  in  the  preceding  sec- 
tions deal  with  the  computation  of  satellite  sub- 
points  relative  to  the  time  and  longitude  of  the 
ascending  node.  Application  of  these  concepts 
to  advance  preparation  of  orbit  schedules  re- 
quires, therefore,  prediction  of  the  times  and 
longitudes  of  ascending  nodes  for  future  orbits. 
The  APT  Predict  messages  give  orbital  period 
to  the  nearest  second  and  longitudinal  displace- 
ment per  orbit  to  the  nearest  hundredth  of  a 
degree.  This  precision  is  adequate  to  extrapolate 
a  few  orbits  ahead.  Greater  precision  is  needed 
for  satisfactory  extrapolation  beyond  a  hundred 
orbits  and  can  be  obtained  from  NESS  or  can  be 
gained  by  empirical  adjustments. 

The  orbital  elements  necessary  for  computa- 
tion of  subpoint  locations  for  extrapolated  orbits, 
also  obtainable  from  NESS,  include  the  length 
of  the  semimajor  axis,  eccentricity,  orbital  in- 
clination, argument  of  perigee  at  a  known  time 
or  reference  orbit,  and  the  rate  of  change  of 
perigee.  The  latter  is  equal  to  the  rate  of  change 
of  a,„  the  value  of  the  angular  displacement  a 
at  the  ascending  node,  which  can  be  computed 
from  the  following  expression  (Runcorn,  1967): 


(/  a„/dt  =  4.98(3M)3.^    (1  ^e2) 
(degrees  per  day). 


^  (5  cos^  i—l) 


where  e  is  the  eccentricity  of  the  orbital  ellipse, 
i  is  the  inclination  of  the  orbital  plane  to 

the  earth's  equatorial  plane, 
a   is  the   earth's   equatorial   radius,   and 
A  is  the  semimajor  axis  of  the  ellipse,  as 

previously  defined. 


For  the  satellite  ESSA-8,  the  value  of  /  is  78.3° , 
A  =  7815.37  km,  and  e  =  0.00323.  With  these 
values  and  taking  6378  km  for  a,  we  compute 


an 


^T—  =    —  1.943  deg  per  day. 

Given  the  nodal  period  of  ESSA-8  as   115.703 
min  we  find,  alternately,  that 


d  «( 
dt 


0.15476  deg  per  orbit. 


The  minus  sign  indicates  that  the  change  is  op- 
posite to  the  angular  motion  of  the  satellite. 

Having  determined  a„  for  a  selected  orbit  and 
starting  at  the  ascending  node  where  a  =a,„ 
the  radial  distance  R  can  be  computed  from  Equa- 
tion (10)  for  use  in  the  numerical  integration  of 
Equation  (9),  to  yield  the  angular  displacement 
(  a  —  a„)  for  successive  intervals  of  desired 
length  after  the  time  of  ascending  node.  The 
values  of  (^  and  A  can  then  be  obtained  with  Equa- 
tions (11)  and  (12).  These  are  the  spherical  co- 
ordinates of  the  point  of  the  orbit  relative  to  a 
nonrotating  frame  of  reference  with  origin  at  the 
earth's  center. 

In  order  to  plot  the  subpoint  track  on  a  chart, 
the  spherical  coordinates  of  points  on  the  orbit 
need  to  be  transformed  to  the  geographic  (geo- 
detic) latitudes  and  longitudes  of  the  correspond- 
ing subpoints.  The  latter  are  the  locations  on  the 
earth's  surface  where  the  local  vertical  passes 
through  the  satellite.  The  angle  which  local  verti- 
cal makes  with  the  equatorial  plane  defines  the 
geodetic  latitude  of  the  subpoint  (Bowditch, 
1958).  Owing  to  the  earth's  ellipsoidal  shape,  a 
line  from  the  satellite  to  the  center  of  the  earth 
intersects  the  surface  at  a  point  north  of  the 
satellite  subpoint.  The  geodetic  latitude  at  this 
point  of  intersection  is  a  close  approximation  to 
that  of  the  subpoint  and  can  be  computed  from 
the  following  expression: 

geodetic  latitude  =  arc  ctn  [(1  ~  q)  ctn  4)] 

where  q  is  the  ratio  of  centripetal  acceleration 
at  the  equator  to  gravity,  as  defined  earlier.  This 
expression  is  derived  from  the  vector  difference 
between  the  true  force  of  gravitation  on  a  unit 
mass  at  a  point  on  the  earth  and  the  apparent 


force  of  gravity  at  that  point,  as  a  result  of  the 
earth's  rotation.  Taking 

q  =  3.468    X  10  "' 

the  maximum  correction  for  geodetic  latitude 
computed  with  the  above  expression  is  about  0.1 
degree  at  latitude  45  degrees. 

The  longitude  of  the  satellite  subpoint  can  be 
obtained  by  adding  the  longitude  of  the  ascend- 
ing node  to  the  value  of  A  computed  with  Equa- 
tion (12).  Adjustment  for  the  earth's  rotation  is 
made  by  adding  to  this  sum  the  product  of  the 
rate  of  rotation  (0.25  degree  of  longitude  per  min- 
ute) and  the  elapsed  time  from  the  ascending 
node  to  the  point  of  the  orbit  for  which  the  sub- 
point  location  is  desired.  A  correction  must  be 
made  for  the  precession  of  the  orbital  plane 
which,  in  a  sun-synchronous  orbit,  amounts  to 
360  degrees  in  365  days  and  is  subtractive. 

APPLICATION 

The  foregoing  procedures  have  been  imple- 
mented in  a  computer  program  for  preparing  or- 
bit schedules  at  the  Southwest  Fisheries  Center, 
La  Jolla  Laboratory,  National  Marine  Fisheries 
Service.  The  program  includes  a  section  for  each 
active  satellite  in  which  appropriate  values  are 
assigned  for  semimajor  axis,  eccentricity,  inclina- 
tion, anomalistic  period,  argument  of  perigee  at 
a  known  reference  orbit  and  the  rate  of  change 
of  the  argument  of  perigee.  All  of  these  quanti- 
ties are  taken  directly  from  information  sheets 
furnished  by  NESS. 

In  running  the  program,  a  data  card  is  submit- 
ted for  each  schedule  to  be  printed.  On  it  are 
punched  the  name  of  the  satellite  and  the  orbit 
number,  day,  hour,  minute,  second,  and  longi- 
tude of  the  ascending  node  for  a  particular  orbit 
as  read  from  a  recent  or  current  APT  Predict 
message.  The  card  is  also  punched  with  the  stand- 
ard (or  daylight)  time  zone  at  the  station  for  con- 
version from  universal  to  local  time,  and  a  speci- 
fication of  the  number  of  minutes  after  ascending 
node  to  which  computed  orbit  reference  times 
will  apply.  The  latter  tells  the  APT  ground  station 
operator  when  to  expect  to  begin  receiving  trans- 
mission. The  number  specified  indicates  to  the 
program  whether  orbits  selected  for  subpoint 
computation  are  to  be  determined  by  south-to- 
north  or  north-to-south  traverses  over  the  station. 
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TABLE  l.—Ver if i cation  of  computed  subpoint  locations  for  ESSA-8.     Computed  values  minus  APT  Predict  values 

are  given  in  degi-ees  latitude  or  longitude. 


Min 

AAN* 

Approximate 
latitude 

Orbit  15239 
(4/12/72) 

Orbit  15641 
(5/14/72) 

Orbit  16030 
(6/14/72) 

Orbit  16444 
(7/17/72) 

Oi-bit  16846 
(8/18/72) 

Orbit  17235 
(9/18/72) 

Orbit  17637 
(10/20/72) 

Lat 

Long 

Lat 

Long 

Lat 

Long 

Lat 

Long 

Lat 

Long 

Lat 

Long 

Lat 

Long 

34 

70° 

.05 

.00 

.04 

-.01 

.06 

.00 

-.01 

-.02 

.04 

-.02 

.00 

.06 

.00 

-.03 

36 

64° 

.03 

-.03 

.06 

.02 

.05 

-.01 

.01 

-.01 

.03 

.00 

-.02 

-.02 

.09 

.00 

38 

59° 

.00 

.03 

.06 

-.04 

.01 

-.01 

-.01 

.07 

-.03 

-.02 

.00 

-.01 

.06 

-.05 

40 

53° 

.01 

.06 

.01 

-.01 

-.01 

.01 

-.02 

.00 

.01 

.00 

.01 

.01 

.07 

-.03 

42 

47° 

.01 

.07 

.05 

-.02 

-.04 

.07 

.05 

.03 

.04 

.02 

.00 

.06 

.08 

-.03 

44 

41° 

.05 

-.03 

.03 

-.01 

-.03 

.01 

.04 

.02 

-.01 

.01 

.03 

.00 

.02 

-.03 

46 

35° 

.04 

-.02 

.06 

-.01 

.03 

.05 

-.01 

-.01 

.00 

-.02 

.01 

.04 

.02 

-.02 

48 

29° 

.01 

.02 

.07 

.03 

.06 

.00 

-.01 

.07 

.06 

.06 

.05 

.00 

.10 

.01 

50 

22° 

.05 

-.03 

.06 

-.02 

.06 

.07 

.07 

.05 

-.01 

.04 

.08 

.06 

.05 

.06 

52 

16° 

.09 

-.01 

.03 

.00 

.05 

.00 

.02 

-.01 

.01 

-.02 

.09 

.00 

.09 

.02 

54 

10° 

.12 

.03 

.10 

.03 

.02 

.04 

.06 

.04 

.01 

.04 

.09 

.04 

.03 

.02 

56 

4° 

.04 

.02 

.06 

.01 

.08 

.02 

.08 

.04 

.00 

.04 

.09 

.04 

.06 

.01 

•"Minutes  after  ascending  node. 


A  sample  schedule,  computed  and  printed  for 
ESSA-8,  is  shown  in  Figure  2.  ESSA-8  transmits 
only  during  the  daytime,  when  it  is  following  a 
north-to-south  track.  Accordingly,  the  reference 
time  appearing  to  the  left  of  each  set  of  subpoint 
tabulations  is  36  min  later  than  the  time  of  as- 
cending node  for  the  orbit  indicated.  The  num- 
bers labeled  "INITIAL  DATA"  were  read  from 
the  data  card  and  are  printed  for  verification. 
The  subpoint  locations  are  printed  to  the  nearest 
hundredth  of  a  degree  of  latitude  or  longitude, 
for  each  minute  from  37  to  56  min  after  ascending 
node,  which  covers  the  daytime  track  from  about 
lat.  60'  N  to  lat.  4"  N.  The  subpoint  location  cor- 
responding to  a  picture  transmitted  by  the  satel- 
lite can  be  found  by  interpolation,  based  on  de- 
termination of  the  time  of  the  picture  relative  to 
the  time  of  ascending  node.  The  .schedule  includes 
only  those  orbits  within  range  of  the  APT  ground 
station  at  the  La  Jolla  Laboratory. 

A  comparison  of  computed  subpoint  locations 
with  those  given  in  APT  Predict  messages  for 
seven  orbits  of  ESSA-8,  completed  between  April 
and  October  1972,  is  presented  in  Table  L  Sub- 


point  locations  are  given  to  the  nearest  tenth  of 
a  degree  in  APT  Predict  messages,  while  com- 
puted values  were  carried  to  the  nearest  hun- 
dredth of  a  degree.  Differences  between  com- 
puted values  and  corresponding  APT  Predict 
values  averaged  less  than  0.06  degree  of  latitude 
or  longitude.  Since  ESSA-8  completes  an  orbit 
in  a  little  less  than  2  hr,  the  motion  of  the  satel- 
lite subpoint  is  about  3  nautical  miles  per  second. 
An  error  of  0.05  degree  latitude  (equal  to  3  nauti- 
cal miles)  in  subpoint  location  is  equivalent, 
therefore,  to  an  error  of  about  1  sec  in  determin- 
ing the  time  of  an  observation  from  the  satellite. 
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